home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / clantp.z / clantp
Encoding:
Text File  |  2002-10-03  |  4.3 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAANNNNTTTTPPPP((((3333SSSS))))                                                          CCCCLLLLAAAANNNNTTTTPPPP((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLANTP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a triangular
  11.      matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      REAL FUNCTION CLANTP( NORM, UPLO, DIAG, N, AP, WORK )
  15.  
  16.          CHARACTER DIAG, NORM, UPLO
  17.  
  18.          INTEGER   N
  19.  
  20.          REAL      WORK( * )
  21.  
  22.          COMPLEX   AP( * )
  23.  
  24. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  25.      These routines are part of the SCSL Scientific Library and can be loaded
  26.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  27.      directs the linker to use the multi-processor version of the library.
  28.  
  29.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  30.      4 bytes (32 bits). Another version of SCSL is available in which integers
  31.      are 8 bytes (64 bits).  This version allows the user access to larger
  32.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  33.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  34.      only one of the two versions; 4-byte integer and 8-byte integer library
  35.      calls cannot be mixed.
  36.  
  37. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  38.      CLANTP returns the value of the one norm, or the Frobenius norm, or the
  39.      infinity norm, or the element of largest absolute value of a triangular
  40.      matrix A, supplied in packed form.
  41.  
  42. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  43.      CLANTP returns the value
  44.  
  45.         CLANTP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  46.                  (
  47.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  48.                  (
  49.                  ( normI(A),         NORM = 'I' or 'i'
  50.                  (
  51.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  52.  
  53.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  54.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  55.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  56.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAANNNNTTTTPPPP((((3333SSSS))))                                                          CCCCLLLLAAAANNNNTTTTPPPP((((3333SSSS))))
  71.  
  72.  
  73.  
  74. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  75.      NORM    (input) CHARACTER*1
  76.              Specifies the value to be returned in CLANTP as described above.
  77.  
  78.      UPLO    (input) CHARACTER*1
  79.              Specifies whether the matrix A is upper or lower triangular.  =
  80.              'U':  Upper triangular
  81.              = 'L':  Lower triangular
  82.  
  83.      DIAG    (input) CHARACTER*1
  84.              Specifies whether or not the matrix A is unit triangular.  = 'N':
  85.              Non-unit triangular
  86.              = 'U':  Unit triangular
  87.  
  88.      N       (input) INTEGER
  89.              The order of the matrix A.  N >= 0.  When N = 0, CLANTP is set to
  90.              zero.
  91.  
  92.      AP      (input) COMPLEX array, dimension (N*(N+1)/2)
  93.              The upper or lower triangular matrix A, packed columnwise in a
  94.              linear array.  The j-th column of A is stored in the array AP as
  95.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  96.              if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.  Note
  97.              that when DIAG = 'U', the elements of the array AP corresponding
  98.              to the diagonal elements of the matrix A are not referenced, but
  99.              are assumed to be one.
  100.  
  101.      WORK    (workspace) REAL array, dimension (LWORK),
  102.              where LWORK >= N when NORM = 'I'; otherwise, WORK is not
  103.              referenced.
  104.  
  105. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  106.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  107.  
  108.      This man page is available only online.
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.